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Abstract. In a follow-up study to a previous analysis of the IRAS 1.2Jy catalogue, we quantify the morphological 
fluctuations in the PSCz survey. We use a variety of measures, among them the family of scalar Minkowski 
functionals. We confirm the existence of significant fluctuations that are discernible in volume-limited samples 
out to 200/i~^Mpc. In contrast to earlier findings, comparisons with cosmological A^'-body simulations reveal that 
the observed fluctuations roughly agree with the cosmic variance found in corresponding mock samples. While 
two-point measures, e.g. the variance of count-in-cells, fluctuate only mildly, the fluctuations in the morphology 
on large scales indicate the presence of coherent structures that are at least as large as the sample. 
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1. Introduction 

Nowadays, one of the primary goals of cosmology is the 
verification of the plethora of models that have emerged 
over the years. To that end, experiments are being con- 
ducted and planned that will ultimately yield deep and 
wide galaxy surveys and high-resolution maps of the 
Cosmic Microwave Background. In parallel to the prepa- 
ration of these large datasets, statistical methods that are 
both easy to interpret and efficient to implement are being 
developed. 

A few years ago, the Minkowski functionals (MFs) 
were introduced into cosmology by Mecke et al. (199^ ). 
The MFs can be calculated from contemporary datasets 
rath er efficiently, an d consequently numerous applications 
(see Kerscher 2000 an d references therein, v ector-valued 
MFs are described in Bcisbart ct al. 2001 ) have estab- 
lished the MFs as a suitable tool for quantifying the mor- 
phological properties of the large-scale structure in the 
Universe . Apart from their practical advantages, MFs 
provide a mathematically well-founded systematic frame- 
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work for the study of statistical morphology and contain 
other higher-order statistics that were employed previ- 
ously. 

The present work extends the scope of an earlier pa- 
per on the morphological fluc tuations in the IRA S 1.2Jy 
catalogue of infrared galaxies ( Kerscher et al. 1998 ), by us- 
ing its much deeper successor, the recently compiled PSCz 
survey (Saunders et al. 2000). However, we do not simply 
repeat the previous study with the enlarged dataset, but 
devise additional sensible tests in order to assess the qual- 
ity of the data and to compare the observational data to 
TV-body simulations. Kerscher et al. (1998D detected, quite 
unexpectedly at that time, significant fluctuations in the 
morphological properties between the northern and south- 
ern hemispheres of the IRAS 1.2Jy catalogue in volume- 
limited samples of up to 200/i~^Mpc depth. While these 
findings were shown to be significant and not due to defi- 
ciencies of the data, simulations could not reproduce them, 
even though they were specifically designed for that pur- 
pose. Meanwhile, simulations have improved, mainly with 
respect to the size of the simulation box. We shall demon- 
strate that the large fluctuations seen in the IRAS 1.2Jy 
catalogue can still be found in the PSCz survey. However, 
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Table 1. Minkowski functionals in three-dimensional 
space expressed in terms of more familiar geometric quan- 
tities. 



Table 2. Volume-limited samples from the PSCz survey. 
N denotes the number of galaxies. 
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state-of-the-art simulations of Cold Dark Matter (CDM) 
scenarios are now capable of predicting sufficient cos- 
mic variance. Another analysis of the PSCz sample with 
Minkowski functionals by Basilakos et al. (2001) focused 
on the shape of super-clusters. 

This article is organised as follows: Sect. ^ briefly ex- 
plains how Minkowski functionals are calculated from the 
PSCz dataset and summarises our findings in a qualitative 
way. It also contains a discussion of possible systematic ef- 
fects and results obtained with other statistical methods. 
Sect. ^ compares the observational data with an analyt- 
ical model and a set of mock catalogues. We summarise 
and provide an outlook in Sect. |[ Finally, a few necessary 
mathematical results are given in the Appendix. 



(i-e-"; 

e~'', 

e-" (1-1^77), 



■q- 



(2) 



Obviously, the Minkowski functionals A* = lj2,3, 

contain a damping factor e^^'. A similar exponential de- 
cay is also found for more general cluster processes. To 
facilitate the comparison between the data and different 
models on large scales, we reduce the functionals by this 
factor and consider the quantities 



MA) 

$P(r) 



(3) 



throughout the following analysis. 

Effective methods to calculate the Minkowski functionals 



2. MepsuHng the morphology of the PSCz survey 



for empirical data points are described in Mccke et al 



(1994 ). To correct for boundary effects without apply- 



2.1. Minkowski functionals 

In order to quantify the morphology within the PSCz sur- 
vey (Saunders ct al. 2000) with Minkowski functionals. 



ing unnecessary statistical assumptions we proceed as 



Schmalzing et al. (1996) (see also Mccke and Wagner 1991 
Kerscher et al. 19971). 



we interpret the redshift space positions {xj^^ of the 2,2. Selection of the samples 



N galaxies in the sample as a realisation of a station- 
ary point process. Adopting the Boolean grain method 
of Mecke et al. (1994), we decorate each point X; with 
a sphere Br(xi) of radius r and consider the union set 



A detailed de scription of the PSCz ga laxy catalogue may 



be found in Saunders ct al. (2000| ). For our analysis 



Ar = [JiLoBr{xi). Hadwiger (1957) has shown that in 
three-dimensional space the four Minkowski functionals 
AIf^{Ar), with /i = 0, 1, 2, 3 provide a complete morpholog- 
ical characterisation of the body Ar . Table |l| summarises 
the interpretations of these functionals in terms of well- 
known geometrical and topological quantities. 

Reduced, dimensionless Minkowski functionals $^(^1^) 
can be constructed by normalising with the Minkowski 
functionals Mf^{Br) of a single sphere. 



(1) 



The mean number density p can be estimated in a sam- 
ple by iV/|fl|, where |f2| is the volume of the sample. 
For a Poisson process these functionals can be calculated 
analytically ( Mecke and Wagner 1991 ) as functions of 
77 = Mf){Br)p, the expected number of galaxies in a sphere 



we extract volume-limited samples (Table ^) from the 
PSCz survey within the standard masked area. The ve- 
locities of the galaxies in the local group frame {v = 
Vhci + sin(/) cos(6) 300km/s) are converted into luminosity 
distances assuming qq = 0.5 and A = 0. 

In order to simplify the boundary corrections, we ap- 
proximate the sample geometry by a spherical cap with 
galactic latitude 6 > 5° for the northern sample and with 
b < —5° for the southern part. Since by doing so, we ne- 
glect some regions which were excluded due to galactic 
absorption or confusion in the IRAS PSC maps, we filled 
these empty regions with random points with the same 
number density as in the fully sampled region. Differences 
in the Minkowski functionals between the filled and un- 
filled samples are barely visible. 

Another issue to be discussed is the flux limit of the 
PSCz survey. Tentatively, we evaluated the Minkowski 
functionals for a series of volume-limited samples with 
100/i~^Mpc depth , but with varying limiting flux. Only 
for limiting fluxes above 0.8Jy the Minkowski functionals 
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stabilise. Likewise, Tadros et al. (1999) find stable results 
only for a flux limit larger than 0.75Jy. Therefore, we shall 
adopt a flux limit 0.8Jy for all of our analyses. 

This kind of flux dependence does not show up in 
two-point measures. The two-point correlation function 
^2, e-g-, does not change significantly when the flux cut 
is varied. Even the mark correlation functions (Beisbart 
and K( Tscher 2000| ; [Szapudi et al. 2000 ) do not show any 
luminosity-dependent clustering at the two-point level. 
Further tests on selection effects will be discussed in 
Sect. |2l. 



2.3. Fluctuating morphology 

Fig. [| displays the values of the reduced Minkowski func- 
tionals 0^ as functions of the rj for the southern and 
northern parts of the PSCz survey at a limiting depth 
of 100/i^^Mpc. For comparison, we also show the expec- 
tation values for a Poisson process with the same number 
density. 

In both parts of the sample the galaxy clustering is no- 
ticeably stronger than in the case of randomly distributed 
points. Moreover, the northern and southern parts differ 
significantly in their morphological features, the north- 
ern part being less lumpy. The most conspicuous features 
are the enhanced surface area in the southern part 
on scales larger than 10/i~^Mpc, and the decrease of the 
Euler characteristic 03, which sets in at 12/i~^Mpc. On 
scales above 10ft.~^Mpc, the integrated mean curvature 
(j)2 and the Euler characteristic 03 are negative, indicating 
an interconnected system of tunnels. In the southern part 
03 remains negative out to large scales, and the large- 
scale structure is dominated by interconnected tunnels 
contrary to the northern part, where completely enclosed 
voids yield positive contributions to the Euler character- 
istic. 



Th o! morphological fluctuations observed in the PSCz 



survey ' appear comparable to thoac found in the IRAS 
1.2Jy catalogue (Fisher et al. 1995; Kerscher et al. 199^ ). 
To verify this, we extract volume-limited samples from 
the PSCz survey with a limiting flux of 1.2Jy. For this 
comparison we use the Euclidean redshift space distance 
as in our original analysis of the IRAS 1.2Jy catalogue. In 
the southern part, we find no difference in the Minkowski 
functionals. The northern part of this PSCz sample, how- 
ever, contains less galaxies than the corresponding IRAS 
1.2Jy subsample, probably due to the more restrictive se- 
lection criteria imposed on the PSCz survey. But still, the 
MFs of the randomly sub-sampled northern part of the 
IRAS 1.2Jy nearly overlap with the PSCz sample. Hence 
we are confident that the MFs pick up the underlying 
large-scale structure, which should of course be the same 
in both catalogues. 



2.4. Selection effects and deeper samples 

To strengthen our claims regarding the morphological fluc- 
tuations, we have to exclude some possible sources of error. 
To begin with, an analysis of the 100/i~^Mpc sample with 
a more conservative cut at |6| > 10° instead of 5° leads to 
nearly identical results for the Minkowski functionals. 

To check whether our results are influenced by some 
peculiar alignment of points, we considered a subsam- 
ple with only 90% of the galaxies randomly chosen from 
the volume-limited sample. Repeating this jackknife pro- 
cedure, we calculate the mean and the error of the 
Minkowski functionals. Still the same prominent fluctu- 
ations between north and south can be seen. These errors 
are smaller than the fluctuations in the random point set. 
We want to emphasise that this jackknife error can only 
serve as an internal consistency check; it does not give 
a reliable estimate of the fluctuations in the underlying 
galaxy distribution ( ^nethlage 199E ). 

Moreover, we select subsamples with "warm" and 
"cold" galaxies as determined from the flux ratio < 2 

and > 2, respectively, with /loo the flux at lOO^m 
and /eo the flux at 60/xm. In both samples we find fluc- 
tuations between north and south, comparable to those 
in the combined sample. Finally, the strength of the fluc- 
tuations hardly diminishes if we consider the measured 
velocities without a heliocentric correction. 

To see whether these fluctuations persist on even 
larger scales, we investigate volume-limited samples with 
200/i-iMpc depth (Fig. |). On smaU scales, both in the 
northern and the southern part, the MFs indicate clus- 
tering. Again and in particular on large scales, the MFs 
of the southern part differ from the MFs of the northern 
part, showing approximately the same pattern as in the 
volume-limited sample with 100ft.~"'^Mpc. The fluctuations 
are now only slightly larger than the one-sigma fluctua- 
tions observed in a Poisson process with the same num- 



ber density. As already demonstrated by Kerscher et al 



(1998), the morphology of such a sparsely sampled cat- 



alogue tends to resemble a pure random sample as dilu- 
tion increases (see also Sect. p?^ ). The Minkowski func- 
tionals in the volume-limited samples with 150/i~"'^Mpc 
depth give comparable results, showing fluctuations of a 
strength between the fluctuations of the 100/i~^Mpc and 
the 200ft.^^Mpc sample. 



In Kerscher et al. (1996) we could show that the mor- 
phological discrepancies in the 1.2Jy galaxy catalogue are 
not a special north-south anisotropy. We proceed similar 
by cutting both the northern and the southern part along 
the yz-plane (in galactic coordinates) into two pieces with 
the same volume. The Minkowski functionals of the four 
samples are shown in Fig. ^ conflrming that the fluctua- 
tions are not a peculiar north-south anisotropy, but actu- 
ally reflect the generic cosmic variance. 
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Fig. 1. The MFs of volume-limited samples from the PSCz survey with 100ft~"'^Mpc depth. The southern (solid) and 
northern parts (dashed) are plotted separately, and for comparison the one-sigma fluctuations (shaded regions) for a 
Poisson process with the same number density are also shown. 




7] 7] 7] 




Fig. 3. The MFs of volume-limited samples with 100/i~^Mpc depth from the PSCz survey. Four different parts of the 
sample are considered separately: southern left (solid), southern right (dotted), northern left (short dash), northern 
right (long dashed). As before, the area shows a Poisson process. 
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and the ratio suggested by van Liesliout and Baddeley 



J(r) = 



1 - G(r) 



1 - F{r) 

For a Poisson distribution one finds 



G{r) = F{r) = 1 - exp ( y^-^ 



(5) 



(6) 



and consequently J(r) = 1. We use the minus-estimators, 
as discussed by Kerscher et al. (1998), to determine F 



10 

r [Mpc/h] 

Fig. 4. Fluctuation of count-in-cells (r) of the volume- 
limited sample with lOO/i^^Mpc depth from the PSCz sur- 
vey. Plotting styles are the same as in Fig. |l|. 

2.5. Other measures 

Both from a physical and a methodological point of view, 
it is interesting to compare the performance of the MFs 
with other statistics. For instance, the well-known (T^(r) 
measures the fluctuations of the number of galaxies in a 
sphere with radius r in excess of a Poisson process. tT^(r) 
is a pure second-order measure. To calculate cr^(r) for 
the PSCz samples, we estimate the correlation integral 
C(r), the mean number of points in a sphere with radius 
r centred on one galaxy, and use the relation 

Cir) 



- 1. 



(4) 



p 47r/3 r-^ 

Our results do not depend on the estimator for C(r) 
(Kerscher 1999). Fig. ^ illustrates that on the two-point 
level there are no significant fluctuations between north 
and south. 

Hence, the fluctuations observed with the MFs are 
a genuine effect of higher-order correlations, illustrating 
that low-order statistics such as the number density or 
the two-point correlation function may miss global fea- 
tures of the large-scale structure. Such higher-order cor- 
relations may originate from coherent elements like fila- 
ments or walls with an extent comparable to or larger 
than the size of the sample, in our case lOO/i^-'^Mpc. This 
conjecture is also supported by visual inspectio n of the 
large-scale structure (see e.g. Huchra et al. 199(i| ). 

Other methods to characterise the spatial distribution 
of points are 

— the spherical contact distribution F{r), i.e. the distri- 
bution function of the distances r of galaxies from an 
arbitrary point in space, 

— the nearest neighbour distance distribution G{r) de- 
fined as the distribution function of distance r of a 
galaxy to the nearest other galaxy. 



and G from the galaxy distribution. The spherical con- 
tact distribution equals the volume density of the first 
Minkowski functional and one minus the void probability: 
F{r)^Mo{Ar)/\n\ - l-Po(Br). 

In Fig. 1^ we show the results for the spherical con- 
tact distribution F{r), the nearest neighbour distribution 
G{r), and the J(r) function. The fluctuations within the 
PSCz survey in comparison to a Poisson process are most 
prominent in the spherical contact distribution F{r), and 
slightly less pronounced in the J (r) -function. Virtually no 
fluctuations are detected with the nearest neighbour dis- 
tance distribution G{r). The reason is that G{r) focuses on 
the small-scale features, which do not fluctuate strongly 
on the scales probed by the PSCz survey, while F{r) and 
J(r) are capable of tracing the large-scale geometry and 
topology of the galaxy distribution. 

2.6. From small to large scales 

In the preceding sections we saw that the morphology on 
small scales up to 10/i^^Mpc does not show the strong 
fluctuations we observed on larger scales. To understand 
this behaviour in more detail we expand the Minkowski 
functionals ^^{Ar) as a series in 77 = Mo{Br)p around 
zero. Based on the expansion in terms of n-point densities 
( iMecke et al. 1994 Eq. (25)) we obtain 



0p(A) 



'2r3 



1 



The integral together with the ratio in front yields a di- 
mensionless geometric factor depending on the correla- 
tion function ^2(5)- The higher-order terms proportional 
to 77""^ include integrals with intersections of n spheres 
weighted by the n-point densities. For small 77 only the 
two-point correlation function is important, but the MFs 
become increasingly more sensitive to higher-order corre- 
lations with larger rj. A small rj may be obtained either 
from a low number density p or a small radius r of the 
spheres. 

The two-point correlation function of the PSCz survey 
shows only tiny fluctuations between north and south as 
can be deduced from Fig. ^. On small scales with 77 ^ 1 
the numerical values of the Minkowski functionals may be 
reproduced using Eq. (Q) with the observed two-point cor- 
relation function and mean density. This explains why for 
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Fig. 5. The panels show the F{r), G{r), and J{r) function of volume-Umited samples from the PSCz survey with 
100/i~^Mpc depth. Plotting styles are the same as in Fig. |^. 



small radii the morphological fluctuations are negligible. 
The expression also explains why the MFs of dilute 
samples tend to look more like those of a Poisson process. 



only should be considered as a first rough estimate. Also 
the zero variance determined from the PSCz for some radii 
is an artifact of using only two samples. 



3. Comparing data and models 

3.1. Second-order moments of Minkowski functional 

To compare the fluctuations of the observed Minkowski 
functionals with different models we will use the dimen- 



sionless second-order moments (see also Eq. (A.2)) 



iv) 



MoiBr 



(( 



(8) 



where \n\ is the sample volume. The m^^ are estimated 
either on the base of samples from simulations or from the 
real data using the standard estimator for a variance. 



MojBr) 1 



1 ^ 

i=l 



M,, 



(9) 



where is the number of samples, = J2i=i -^/j 

is the sample mean, and the M;} are the MFs of the 
ith sample. From the PSCz catalogue we take two sam- 
ples, the southern and the northern part with a depth of 
lOO/i-^Mpc each. 

In Appendix A the fluctuations for a Poisson and a 
binomial process are given analytically. The comparison 
with the estimated fluctuations from the PSCz survey in 
Fig. ^ indicates that the observed dimensionless fluctua- 
tions are significantly larger than the for a Poisson 
process with the same number density. Qualitatively, the 
fluctuations from the PSCz survey show similar features 
as for a Poisson process, i.e. the functional form and the 
relative strengths of the fluctuations are comparable. Only 
the surface area is qualitatively different on large scales. 
Of course, an estimate of fluctuations from two samples 



3.2. Comparison with CDM simulations 

More realistic models to compare our data to are based on 
A^-body simulations. We use mock galaxy catalogues con- 
structed from the VIRGO consortium's Hubble volume 
simulations. The mock catalogues were kindly provided 
to us by Shaun Cole. From these mock catalogues we ex- 
tract volume-limited samples with a depth of 100/i~^Mpc 
similar to the observed galaxy samples. The CDM A^- 
body simulations were performed by the VIRGO con- 
sortium and comprise a ACDM (il.,„ = 0.3, JIa — 0.7, 
(78 = 0.9) and a rCDM {VL^ = 1, f^A = 0, cts = 0.6) simu- 
lation, conducted in a box with a side length of 3/j-~^Gpc 



and 2h ^Gpc, respectively (for details see Evrard et al 
[2001 and their web-sit(^. Mock catalogues were gener 



ated with methods similar to those used by Cole et al 



(1998): the galaxies were extracted from the distribution 



of the dark matter particles employing a biasing scheme; 
luminosities have been assigned according to the observed 
luminosity function (for details see Cole et al. 199S and 
their web-site^). The observer position was not chosen at 
random but selected to satisfy criteria similar to those of 
the local group, with respect to the local density, the CMB 



dipole, and the local shear (Baugh et al. 2001). We used 
the first eleven mock catalogues constructed from each 
simulation. 

The Minkowski functionals for the mock samples were 
calculated in redshift space within the same sample geom- 
etry as for the PSCz survey. Considering north and south 
separately we can use A^., = 22 samples per simulation 
to estimate the mean MFs and their variance. In Fig. |^ 
the mean MFs of the mock samples are compared with 



tittp : //www .physics . lsa.iimich.edu/hubble-volume/ 



tittp : / / star-www. dur . ac .uk/~cole/mocks/hubble .html 
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Fig. 6. The dimensionless fluctuations of the Minkowski functionals for a Poisson process are shown in the left plot. 
In the middle plot the estimated variances m^^ from the volume-limited sample of the PSCZ with 100ft.~^Mpc depth, 
and in the right plot the estimated variances from the ACDM mock samples are shown (notice the different scaling). 
The rCDM mock samples show similar fluctuations. 
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T] T] 7] 



Fig. 7. In addition to the mean MFs from the PSCz survey (solid), this plot displays the mean results from the mock 
sample constructed from the ACDM (dotted) and rCDM (dashed) simulation. The shaded area marks the one-sigma 
region determined from the ACDM simulation. The tCDM simulation shows similar fluctuations. 



the mean MFs of the PSCz catalogue. Within the esti- 
mated errors, the large scale morphology of the galaxy dis- 
tribution is marginally reproduced by the ACDM model, 
whereas the rCDM model is too lumpy. 

The main issue is whether these mock catalogues are 
able to reproduce the large morphological fluctuations in 
the PSCz survey. This is of specific interest since it turned 
out that a simulation, designed to match the local distri- 
bution of IRAS galaxies, was not able to reproduce the 



fluctuations in the 1.2Jy galaxy survey ( Kerscher et al 
1998|)~|Similar large-scale fluctuations as observed in the 
galaxy distribution may indeed be found in the mock cat- 
alogues of both models (see Fig. ^). At least, the ampli- 
tude of the fluctuations m^^ determined from real data 
and the mock samples are of the same order of magni- 
tude. The fluctuations in the number density between of 
the volume-limited samples extracted from the mock cat- 



alogues are in full agreement with the cr^ expected for the 
models at the sample size, but are larger than the number 
density fluctuations between the northern and southern 
PSCz samples. These fluctuations in the number density 
seem to cause the increased fluctuations in moo for the 
mock samples. 



4. Summary 

We investigate the morphology of the galaxy distribution 
in the PSCz survey. In order to quantify the geometri- 
cal and topological features, we use the fluctuations of 
count-in-cells as a two-point measure, as well as a number 
of higher-order statistics, namely the nearest neighbour 
distance distribution function, the void probability, the 
J-function, and the family of scalar Minkowski function- 
als. We find appreciable fluctuations between the northern 
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and southern parts of the survey, confirming earlier find- 



ings frcim the IRAS 1.2Jy catalogue ( Kerscher et al. 1998|) 
These fluctuations are discernible in volume-limited sam 
pies of as much as 200/i^^Mpc depth. By inspecting sub- 
samples we show that these fluctuations are not due to a 
specific north-south anisotropy, but are an intrinsic fea- 
ture of the galaxy distribution on large scales. 

The comparison with the fluctuations in mock- 
catalogues shows that iV-body simulations of two CDM 
models in boxes with a side length of 2 and 3/i^^Gpc fol- 
lowed by a suitable biasing procedure are able to account 
for the observed morphological fluctuations in the galaxy 
distribution. This was not the case when we compared the 
fluctuations in the IRAS 1.2Jy with those determined from 
a simulation in a smaller box with 250ft.^^Mpc side length. 
Concerning the mean values of the Minkowski function- 
als, the rCDM simulation is not able to explain the ob- 
served large-scale morphology. Marginally, the ACDM 
model succeeds. 

Speaking of fluctuations, one usually thinks of the fluc- 
tuations in the number of galaxies, the fluctuations of 
count-in-cells cr^, or the power spectrum quantifying the 
amount of fluctuations in a Fourier mode of the density 
contrast. Both measures quantify the galaxy distribution 
at the two-point level, and are therefore unable to de- 
tect coherent large-scale structure elements like walls and 
filaments (sec e.g. Szalay 1997; Bchmalzing et al. 1999; 



A^-body simulations could verify (Melott and Shandarin 
199C; Jenkins et al. 1998| ). To detect these morphological 



Chiang and Coles 2000). Indeed, our investigations of the 
galaxy distribution in the PSCz survey using show no 
significant fiuctuations at the two-point level between the 
northern and southern parts of the galaxy sample. 

Minkowski functionals, however, incorporate correla- 
tions of arbitrary order and quantify the connectivity 
of large-scale structure elements in the galaxy distribu- 
tion. The amplitude of the morphological fiuctuations in 
samples with an extent of 100/i^^Mpc is approximately 
three times larger than expected for a Poisson process. 
The morphological fluctuations remain detectable even 
in a 200/i~^Mpc sample, but since the deeper volume- 
limited samples are sparser and therefore more similar to 
a Poisson process, the results become less signiflcant. The 
structural differences between the northern and southern 
parts of the PSCz sample show up prominently in the 
Minkowski functionals for radii above lO/i^^Mpc. As al- 
ready mentioned, the number density and also the two- 
point correlations show only small fluctuations. This in- 
dicates that a variety of coherent structures with a low 
density contrast and an extent of more than lOO/i-^Mpc, 
perhaps 200/i~-'^Mpc, are shaping the large-scale structure 
of the Universe. In this sense, our analysis quantifles the 
visual impression of features like the great wall in the CfA2 
catalogue (Huchra et al. 1990). 

Structures like walls and filaments were predicted by 
analytical an d numerical work based on the Zel'dovich ap- 
proximation ( Arnol'd et al. 1982 ; Doroshkevich et al. 1996 ) 



and re late d approximations ( Kofman et al. 1992| ; Bond 
et al. l |996|) . They are generic features of the gravitational 
collapse for CDM-like initial conditions, as numerous 



features, especially if they show up only at a low density 
contrast, sensitive methods like the Minkowski functionals 
are needed. They allow us to quantify the geometry and 
topology of the network of walls, filaments, and voids seen 
in the distribution of galaxies. 
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Appendix A: Fluctuations of the Minkowski functionals for a Poisson model 

The simplest point process, and thus the simplest reference model to compare real data to, is a stationary Poisson 
point process lacking any spatial correlations. In this appendix we give the exact expressions for the variances of the 
Minkowski functionals for a Poisson process and a binomial process. These results prove useful both to estimate the 
fluctuations of the morphological properties for finite samples, and to serve as a reference model for a comparison 
with fluctuations in empirical datasets (Fig. ^). Here we only explain the basic notions and list the results for three 
dimensions. A full derivation of the variances for the Minkowski functionals as well as their covariances in any dimension 



can be found in (Mecke 2001) 



Let us first consider a stationary Poisson point process with number density p in a sample Q of volume \Q\ 
embedded into three-dimensional Euclidean space. We are interested in the morphological fluctuations as quantified 
by the covariances of the Minkowski functionals Mf^{Ar) for the Boolean grain model. For convenience, we divide the 
variances by the Minkowski functionals M^{Br) of a single sphere of radius r and express them as functions of 
77 = MQ(Br)p such that all quantities are dimensionless and may be compared quantitatively: 

where (•) denotes the average over the statistical ensemble. 

The Minkowski functionals Mfj.{Ar) of a union set of balls can be represented either as an integral over the volume 
of the set (/i = 0) or as integrals over the {d — A)-dimensional intersections of A spherical boundaries for fi — I, . . . ,d 
and A — 1,... in the latter case local curvatures arise as weightings. Therefore, the variances Ti^j^^ can be 
decomposed into a series of curvature-weighted structure funct ions M^lj^^ ( ri, s) containing the contributions of the 



{d — Ai) -dimensional spherical intersections to Mfj_. for i = 1,2 ( Mecke 2001 ) 



m,,M - E E / ^^Mi'M.'C'y;^)- (A.2) 

For the following listing of the structure functions V denotes the volume of the sphere Br, and V{s) is the normalized 
volume of two overlapping spheres i?f.(xi) and _Br(x2) of radius r and centre given by 

m(xr)Ui?.(x.))^ 3 1 
^ ' V{Br) 2 2' ^ ' 

where s = ^^^2r^^^ denotes the normalized distance of the spheres. Note, that the structure functions arc symmetric, 
M^i^2 g-j = M^^Ai (i^^ vanish, if A^ > /i^ for i = 1 or i = 2. We get 

M°°{ri; s) = 24s2 (^c-'^v(s) _ ^-2^^ ^ (^ 4) 

^2 „Ai + A2~i 

/lfAiA2C . \ /'_iUi + A2„-))y(s) / Af-^1^2 _ 24s2 



The structure amplitudes M^, M^{f]), M^{s) and M^^^l^is) arising in these formulae vanish if A > /i and A^ > fii for 
i = 1 or i = 2, respectively. For the other cases one finds the following: the quantities M^, with 



2Mi{Brf 37r2 



Mi 



M3 



Mo{Br)M2{Br) 16 
(Br) . 

3 _ 6Mi{Brf _ 977^ 

Mo{Br)Ul3{Br) ~ Te" 



- 2 ^ 6Afi(i3,)M2(i?.) ^ (A.5) 

3 M,{Br)Ah{Br) ' 



contain local curvature contributions. The M^(rj) additionally incorporate the probability of finding A points within an 
otherwise empty sphere B^ resulting in an intersection which - according to the additivity of the Minkowski functionals 
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27r 



da;idx2y 1 - ixiX2 + Jl - xfJl - cos^)^ F^[t] 



2 J Mo{Br) 
-s 

1 27r 1 

MI{S) = ^ Jdx JJ-^^^ jjdXidX2A[xi,X2,X,(l>i,(l>2] 



(A.6) 



-s 



x\/ 1 — x\\ 1 — x^ sin(02 — </'i) — xi\/ 1 — x"^ \ 1 — sin 02 + X2 \/ 1 — a;^ a / 1 — sin0i 



Finally, we turn to the M^l^l.i{s). For / = 0, one finds that they split into a product of two functions already 



known from above: 



(A.7) 



For Z 7^ 0, the expressions become more involved: 



fJ.lfJ.2 


lis) 


= 2s, 




lis) 


= 2s 




lis) 


1 

~ 2s 




2is) 


-I 




Mil 


lis) 


= 2s 



Mll.iis 
Mll.,2is 

Mii,,is 

Mii.2is 



2tt 1 

2s J J darF^Ji(-s,ar,<^)]VW^^^+7r^^Vr^^^W?, 



-M^,li,iis)M^l,.^,{s), 

d(l) 



MoiBr) 



1 - (s^ + (1 - s^) COS 0)^ I F^, [t{-s, -s, <I})]F^, [t{-s, -s, 
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\/l — s'^X2 \Jl — x\ sin (j)i — Vl — s'^Xi '^1 — X2 sin (j)2 — 1 — x\^l — X2 sin(02 — (j'l) 



^Mlliis)Mll,.^,{s), 



27r 



= jj Mo{B% J '^''^^ ~ + ~ cos0i)2A[-s, -5, X, <Ai, <^2]F^. [ti-s, -s, h - </>!)] 



(1 — s^)xsin(^i + s-\/l — s^y^i _ ^2 (sin((/)i — ^2) +sin02) 



2s 

2tt 



2Tr 



_^da.,dx2A[a.i,-s,-s, A[a.2,-s,-s,<^2,<^] 

-s 
2 ^ 

<Y\ il- s'^)xi sin + -\/r^-s^Y 1 - a;| (sin(0 - + sin ^j) 



^3^3(3) = 2^(1 - s') 



2tt 

2n2 // d0ld02 



)i# |Sin(0i - (j)2) - sin(/>i +sin(/>2| (A[-s, -s, -s, 1^2]) + 2M (5(xi - X2), 



(A.8) 
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Fig. A.l. The normalized covariances of the Minkowski functionals of a Poisson process (see Fig. |6| for the variances). 



where we have used the functions Fi: 

-F2[i] = arcsini(xi, X2, 0), i^sM 
1 



2t{xi,X2, 



t{xi,X2, 



V2 



(A.9) 



1 - XiX2 



1 — xf\/ 1 — X2 COS ( 



Furthermore, Eq. (A. 8) employs the spherical excess A[xi, X2, x, (f)i, 4)2] given by I'Huilier's formula: 



tan 



A 



tan 



ai + a2 + as 



tan 



ai + a2 — as 



tan 



tan ■ 



-ai +02 + 03 



4 4 4 4 

where the edges of a spherical triangle are given by = 2 arcsin and 

h =t{x2,X,(j)2), t2=t{x,Xi,-(j)i), t3 =t{xi,X2,(f>l - (f>2)- 



(A.IO) 
(A.ll) 



And lastly, the pre-factor of the (5-distribution at the end of Eq. (A. 8) is given by 



2tt 



M- 



W2 



-Idx 



id02 



da;idx2A[xi, X2, 01, 02]" 



(A.12) 



x\/ 1 — xi\/l — X2 sin 



bi) — xi \/l — x'^\J 1 — x\ sin 02 + X2 \/ 1 — x"^ \J 1 — x\ sin 0i 



While the expressions appear somewhat complicated at first glance, they can be easily evaluated by numerical inte- 
gration. 

In Figs. ^ and Al the second-order moments are shown as functions oi rj — Mo{Br)p. One finds an increasing 



behaviour, mi,fj,{vi) = rj + 0(77 ), for low densities 77, and the asymptotic behaviour mi,f^{r] 



3) ^ jy-^+A'-de-^'. The 



Euler characteristic shows by far the largest variances. 

Whereas for a Poisson process (P) the number of points within our sample is fluctuating, a binomial process (B) 
corresponds to a canonical ensemble, where the number of points is fixed. The fiuctuations of the Minkowski functionals 
for a binomial process can be derived from the Poisson case, since one finds that in general (Mecke 2001) 



d7] 



drj 



(A.13) 



where ^^Av) the mean values given by Eq. (Q). Results for a binomial process are shown in Fig. [A.2| . As shown 
in Fig. A. 2, the binomial process exhibits smaller fluctuations compared to the Poisson case depicted in Fig. |A.1| . The 
qualitative features in the fluctuations are similar but significantly damped. Also additional correlations between the 
volume, the surface area, and curvatures show up. 
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Fig. A. 2. The covariances of the Minkowski functionals, now for a binomial process. 



